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I.  INTRODUCTION 


The  use  of  graphitic  materials  for  thermal  protection  system  applica¬ 
tions  involving  high  heat  flux  environments  is  based  upon  the  fact  that  car¬ 
bon  possesses  a  rare  combination  of  very  desirable  properties.  Its  thermal 
properties  make  graphitic  materials  good  heat  sinks  with  a  very  high  melt 
temperature,  while  its  thermochemical  properties  indicate  a  very  large 
heat  of  sublimation.  A  considerable  amount  of  ablation  response  modelling 
for  graphite  and  graphitic  materials  has  been  reported  (Ref.  1). 

Until  recently,  the  maximum  heat  flux  environments  encountered,  in 

ballistic  reentry  vehicle  applications  and  arc  jet  (Ref.  2)  ground  tests, 

2 

have  been  about  10  kW/cm  .  At  the  present  time,  however,  combined  con- 

2 

vective/ radiative  heat  transfer  up  to  50  kW/cm  is  predicted  for  the  Jupiter 
entry  probe  (Ref.  3)  which  will  probably  have  a  graphitic  material  heat 
shield  (Ref.  4)  and  laser /material  interaction  tests  have  exposed  graphite  to 
fluxes  of  100-1000  kW/cm2  (Ref.  5)  and  even  as  high  as  10^-10^  kW/cm2 
(Ref.  6). 


Hurwirz,  H.,  K.M.  Kratsch,  and  J.E.  Rogan,  "Ablation,  "  AGARDograph 
No.  161,  March  1972. 

2 

Lundell,  J.H.,  and  R.  R.  Dickey,  "Ablation  of  ATJ  Graphite  at  High 
Temperatures,"  AIAA  Journal,  11,  February  1973,  pp.  216-222. 

3 

Moss,  J.N.,  E.E.  Anderson,  andC.W.  Bolz,  "Aerothermal  Environment 
for  Jovian  Entry  Probes,  "  AIAA  Paper  76-469,  AIAA  11th  Thermophysics 
Conference,  San  Diego,  Calif.,  July  1976. 

4 

Nunamaker,  R.R.,  quoted  in  Aviation  Week  and  Space  Technology, 

April  3,  1978. 

5 

Lundell,  J.H.,  and  R.  R.  Dickey,  "Radiative  Vaporization  of  Graphite  in 
the  Temperature  Range  of  4000°  to  4500°K,  "  Thermophysics  of  Spacecraft 
and  Outer  Planet  Entry  Probes,  Progress  in  Astronautics  and  Aeronautics, 
S6,  1977. - - - - 

^Stegman,  R.  L.,  J.T.  Schriempf,  and  L.  R.  Hettche,  "Experimental  Studies 
of  Laser -Supported  Absorption  Waves  with  5  ms  Pulses  of  10.  6fi  Radiation,  " 
Journal  of  Applied  Physics,  44(8),  August  1973,  pp.  3675-3681. 


Ablation  response  models  currently  in  use  for  graphite  in  these  very 
high  heat  flux  environments  rely  upon  direct  extension  of  the  classical 
modelling  (Ref.  1)  which  assumes  solid-gas  phase  equilibrium.  Baker  (Ref.  7) 
has  developed  a  nonequilibrium  model  for  graphite  sublimation  in  an  air 
environment  and  concluded  that,  for  relatively  high  heat  flux  at  low  pres¬ 
sures,  nonequilibrium  effects  can  become  important. 

The  primary  purpose  of  this  work  has  been  to  model  the  nonequilibrium 
thermochemical  response  of  graphite  exposed  to  external  radiation  environ- 
ments  extending  to  1000  kW/cm  .  Above  this  flux  level,  laser-supported 
absorption  wave  effects  (Ref.  6)  tend  to  significantly  decouple  the  radiation 
input-material  response  behavior.  The  work  has  proceeded  in  two  stages 
and  is  described  in  this  manner.  In  the  first  stage,  the  Knudsen- Langmuir 
type  model  developed  earlier  (Ref.  7)  was  used.  In  the  second  stage,  an 
irreversible  thermodynamics  phase  change  model  developed  by  Bornhorst 
and  Hatsopoulos  (Ref.  8)  was  adapted  to  the  present  problem  and  extended  to 
include  multiple  species  with  chemical  reactions. 

The  problem  formulation  using  these  two  models  is  given  in  the  next 
section.  Then,  the  model  predictions  are  compared  with  available  experi¬ 
mental  data.  Finally,  some  preliminary  conclusions  and  recommendations 
for  further  work  are  given. 


7 

Baker,  R.  L.,  "Graphite  Sublimation  Chemistry  Nonequilibrium  Effects,  " 
ALAA  Journal,  15(10),  pp.  1091-1097. 

g 

Bornhorst,  W.  J.,  and  G.  N.  Hatsopoulos,  "Analysis  of  a  Phase  Change  by 
the  Methods  of  Irreversible  Thermodynamics,  "  Transactions  of  the  ASME, 
J.  Appl.  Mech.,  34(4),  December  1987,  pp.  840-846. 


-6- 


II.  PROBLEM  FORMULATION 


A.  KNUDSEN- LANGMUIR  MODEL 


The  Knud  sen- Langmuir  equation  is  obtained  from  kinetic  theory  con¬ 
siderations  and  is  given  by  Kennard  (Ref.  9)  as 


m 


(la) 


where  m  is  the  mass  flux  between  phases  per  unit  wall  area,  is  the 
molecular  weight,  T^  is  the  wall  temperature,  pv  is  the  vapor  pressure  of 
the  vaporizing  material  at  the  wall  temperature,  p  is  the  partial  pressure  of 
vaporizing  material  at  the  surface,  and  a  is  the  vaporization  coefficient. 

For  application  to  graphite  ablation,  Eq.  (la)  has  been  assumed  applicable 
to  each  carbon  species  (Refs.  7,  10),  i.  e. 


m.  =  a. 
1  x 


Mi 

27tRT 


(lb) 


The  total  mass  loss  is  then  obtained  by  summing  over  n  species 


m 


(lc) 


9 

Kennard,  E.H.,  Kinetic  Theory  of  Gases,  McGraw-Hill  Book  Company,  Inc., 
New  York,  N.  Y„  1938. 

^Dolton,  J.A.,  R.E.  Mauer,  and  H.E.  Goldstein,  "Thermodynamic  Per¬ 
formance  of  Carbon  in  Hyperthermal  Environments,  "  AIAA  Paper  No. 
68-754,  AIAA  Third  Thermophysics  Conference,  Los  Angeles,  Calif., 

June  1968, 
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Since  the  primary  interest  of  the  present  work  is  high  external  radia¬ 
tion  environments,  diffusive  (boundary  layer)  heat  transfer  effects  will  be 
assumed  to  be  negligible  and  not  retained  in  the  surface  energy  balance.  It  is 
further  assumed  that  the  suriace  mass  addition  rates  are  sufficiently  high  so 
that  only  carbon  species  are  present  at  the  ablating  surface.  With  these 
simplifying  assumptions,  the  surface  energy  balance  becomes  simply 


icond 


=  q  .  -  €(TT  -  m[H  -  H  ] 
^rad  w  1  w  s 


(2a) 


.  .  4 

where  <lcon(j  i-s  the  heat  conducted  into  the  ablating  material,  qra(j  andcaT^ 

are  the  net  heat  radiated  to  and  away  from  the  surface,  and  m[Hw  -  Hg]  is 
the  net  heat  convected  away  from  the  surface  by  the  mass  loss.  Assuming 
the  heat  conducted  into  the  material  is  given  by  the  steady- state  approxima¬ 
tion,  which  is 

q  ,  =  m  c  [T  -  Tj  (2b) 

^cond  pL  w  bJ 


the  surface  energy  balance  equation  becomes,  with  Hg  =  c^Tw 


q  -  €<7T4  -  m[H  -  H  J  =  0 
^rad  w  ^ 


w 


(2c) 


For  the  multiple  species  case,  the  wall  enthalpy  is  given  in  terms 
of  the  partial  specific  enthalpies  h^  and  the  mass  fractions  by 

n 

H  =  V  K.h.  (3) 

w  /  j  11  ' 

i=l 


Also,  the  mass  fraction  IC  is  represented  by 


K. 

l 


V'i 


(4a) 
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and  the  average  molecular  weight  is 


n 


-i 


(4b) 


i  =  l 


i=l 


i=l 


The  Mach  number  of  the  flow  from  the  ablating  surface  is  defined  by 


M 


V 

/  yRT 

V  j 


(5) 


To  determine  the  partial  pressure  of  species  i  at  the  surface,  it  is 
necessary  to  assume  either  frozen  or  equilibrium  chemistry  in  the  gas  phase. 
If  the  chemistry  is  assumed  to  be  frozen,  then 


m. 

i 


=  rnK. 

l 


(6a) 


By  combining  Eqs.  (lb),  (lc),  (4a),  (5)  and  (6a)  we  obtain  for  this  case 


-1 


,  ,  M  , 

Pi  = 

1  +w 
i 

V  J 

(6b) 


if  the  chemistry  is  assumed  to  be  equilibrium  in  the  gas  phase,  then  the 
partial  pressures  pi  are  given  by 

Pi  *  [p^Pp  (6c) 
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The  equilibrium  vapor  pressure  p^  and  the  partial  specific  enthalpies  h^  of 
species  i  at  the  wall  temperature  are  represented  by 


p^  (atm)  = 


Ai 

T{dK)  +  Bi 


(7a) 


h. 

x 


2500 


+  b.T 
1 


dT'  +  d. 


l 


(7b) 


The  values  of  the  constants  in  Eqs.  (7a)  and  (7b)  for  the  three  thermochemi¬ 
cal  data  bases  (Refs,  10-12)  used  in  the  present  study  are  given  by  Baker 
(Ref.  13). 

Equations  (1)  through  (7)  provide  (5n  +  7)  equations  in  the  same  number 
of  unknowns  to  be  solved  simultaneously  to  obtain  mass  loss  rate  m  as  a 
function  of  qracj  and  the  external  pressure  pe<  If  the  Mach  number  calculated 
from  Eq.  (5)  is  supersonic,  then  pfi  is  increased  in  accordance  with  the 
requirement 


Laser/material  interaction  flows  of  this  type  have  recently  been  reported  by 
Covington,  et  al,  (Ref.  14),  who  observed  a  plume  flowfield  away  from  the 
ablating  surface  with  a  barrel  shock  and  a  Mach  disk. 

^"JANAF  Thermochemical  Tables,  "  The  Dow  Chemical  Co.,  Midland, 
Michigan,  1969. 

12 

Palmer,  H.  B.,  and  M.  Shelef,  "Vaporization  of  Carbon,  "  Chemistry  and 
Physics  of  Carbon.  4,  ed.  P.  L.  Walker,  Jr.,  Marcel  Dekker,  Inc.,  New 
York,  N.Y.,  1 968. 

13 

Baker,  R.  L.,  "The  Thermochemical  Response  of  Graphitic  Materials 
Subjected  to  High  External  Radiation  Environments,  "  Report  No.  TOR -0078 
(3701)-2,  The  Aerospace  Corporation,  El  Segundo,  Calif.,  April  1978. 

14 

Covington,  M.A.,  G.  N.  Liu,  and  K.  A.  Lincoln,  "Free-Jet  Expansions  from 
Laser- Vaporized  Planar  Surfaces,  "  AIAA  Paper  No.  76-22,  ALAA  14th 
Aerospace  Sciences  Meeting,  Washington,  D.  C.,  26-28  January  1 976. 


-10- 


B. 


IRREVERSIBLE  THERMODYNAMICS  MODEL- - 
SINGLE  COMPONENT  (C3) 


The  Knudsen- Langmuir  model  discussed  above  considers  the  pressure 
discontinuity  Ap  =  p  -  pv  but  neglects  the  possible  existence  of  a  temperature 
discontinuity  AT  =  T  -  T  .  It  also  neglects  potential  coupling  between  the 
mass  flux  m  and  the  energy  flux  j  ,  where  j  =  H  m.  Irreversible  thermo- 

U  U  W  jJ. 

dynamics  provides  a  methodology  whereby  these  effects  can  be  modelled. 

The  first  application  of  this  methodology  to  phase  change  modelling 
was  apparently  the  work  of  Bornhorst  and  Hatsopoulos  (Ref.  8).  In  the  fol¬ 
lowing  paragraphs,  the  salient  features  of  their  development  are  given.  The 
relationship  of  this  type  of  formulation  to  the  Knudsen- Langmuir  represen¬ 
tation  and  the  groundwork  for  extension  of  the  irreversible  thermodynamics 
model  to  the  multiple  species  chemically  reacting  case  are  then  discussed. 

In  order  to  model  the  discontinuous  behavior  of  intensive  properties  at 
the  subliming  graphite  surface,  the  discontinuous  system  is  envisioned  as 
shown  in  Fig.  1.  Two  equilibrium  regions  I  and  II  are  connected  by  a  region 
of  vanishingly  small  thickness,  III.  The  pressure  and  temperature  in 
regions  I  and  II  are  considered  to  be  held  constant  by  pistons  and  heat  reser¬ 
voirs.  These  quantities  are  highly  nonuniform  in  the  communicating  region 
III. 

The  total  entropy  production  <r  for  the  system  depicted  in  Fig.  1  is 
(Refs.  8,  15) 


tot 


=  -j  x  -  ]  x 
u  u  Jm  m 


(9) 


1 5 

DeGroot,  S.  R.,  and  P.  Mazur,  Non -equilibrium  Thermodynamics,  North 
Holland  Publishing  Co.,  Amsterdam,  19^2. 

jjj 

The  entire  methodology  is  based  upon  the  validity  of  linear  flux-force  rela¬ 
tionships  [see  Eqs.  (11a, b)]. 


RESERVOIR  I 


where  and  jm  are  the  net  interphase  energy  and  mass  transport  rates 
(j  se  m),  and  x^  and  xm  are  the  driving  potentials  or  affinities  for  energy 
and  mass  transfer,  respectively,  and  are  given  by 


•*(*) 

•*(*) 


(10a) 


(10b) 


where//  is  the  chemical  potential. 

The  phenomenological  equations  relate  the  transport  rates  and  driving 
potentials  in  terms  of  the  transport  coefficients  L.  Thus 

-  -  Lu At)  -  LuiA(r)  <Ui 

im  *  -  hAk)  -  LiiA(f )  <UI 

The  chemical  potential  for  a  homogeneous  system  is  related  to  the 
temperature  and  pressure  by  (Ref.  16) 

£  =  h|KT)  +  in  p|  (12; 

It  can  easily  be  shown  that 


Callen,  H.G.,  Thermodynamics,  John  Wiley  &  Sons,  Inc.,  New  York, 


In  Ref.  8  the  transport  coefficients  L  are  evaluated  by  considering 
various  limiting  cases  of  j  or  AT  going  to  zero  and  comparing  the  resultant 
form  of  Eqs.  (11a)  and  (lib)  with  kinetic  theory  expressions  for  that  par¬ 
ticular  limiting  case.  By  combining  the  trarisport  coefficient  expressions 
obtained  in  this  way  with  Eqs.  (11a),  (lib),  (lZd ),  and  the  Onsager  reciprocity 

condition  L  .  *  L.  ,  Bornhorst  and  Hatsopoulos  obtain 
ui  1U 


j 


u 


(h 


1  RT\  .  AT 

2  Jt  I  Jm  -  he  t2 


(13a) 


where 


j  =  L.. 
Jm  u 


pvAT  vAp 
2T2  T 


H  *  \fzdr  pcvtZ 

i  - 1  /  \  /  J 

ii  “  v  \2-a  /  Y  27TRT 


(13b) 


(14a) 

(14b) 


and  v  is  the  specific  volume  which  is  the  reciprocal  of  the  mass  density. 

Of  particular  interest  is  the  limit  of  j  for  AT  =  0. 

(For  AT  =  0)  jm  =  -  <15> 


Comparing  this  expression  with  the  Knudsen- Langmuir  expression  given  by 
Eq.  (la),  we  see  that  they  are  equivalent  except  that  the  vaporization  coeffi¬ 
cient  (Ref.  12)  appears  as  ct  in  Eq.  (la)  and  in  the  term  [(2a/2-a)]  in  Eq.  (15). 


1 


The  reason  for  this  was  first  recognized  by  Schrage  (Ref.  17)  who  demon¬ 
strated,  using  kinetic  theory  arguments,  that  by  accounting  for  the  gas  phase 
bulk  flow  velocity  in  the  velocity  distribution  function  at  the  surface,  one 
obtains,  for  small  j  Eq.  (15)  in  place  of  Eq.  (la).  Physically,  this  indi¬ 
cates  that  the  bulk  flow  velocity  at  the  surface  gives  rise  to  an  effective 
vaporization  coefficient  [  (2ct/2 -Of)].  The  basic  definition  of  01  (Ref.  12), 
however,  has  not  been  changed. 

Two  alternate  forms  of  Eqs.  (13a)  and  (13b)  are  useful  as  an  aid  in 
understanding  the  relationships  of  the  irreversible  thermodynamics  formu¬ 
lation  to  the  Knud  sen- Langmuir  representation  and  of  the  single  component 
model  to  the  multiple  species  chemically  reacting  model.  Thus,  following 
DeGroot  and  Mazur  (Ref.  15)  if  we  require  the  entropy  production,  as  given 
by  Eq.  (9),  to  be  constant  and  define  a  new  energy  flux  j  and  a  volumetric 
flux  jv  by 


jq  “  ^u  "  Hw^i 


Jv  =  VJm 


(16a) 


(16b) 


then  a  symmetric  set  of  equations  entirely  equivalent  to  Eqs.  (13a)  and 
(13b)  is 


i  =  -  A 


-  A 

qv  T 


a  AT  a  Ap 
Jv  =  '  Avq  x2  '  Avv  T 


(17a) 


(17b) 


'Schrage,  R.  W.,  A  Theoretical  Study  of  Interphase  Mass  Transfer,  Columbia 
University  Press,  New  York,  N.  Y.,  1953, 


where 


A 

qq 

-  He  +  7  (P^2  Hi 

(18a) 

A 

qv 

=  A  =  -  jpv2  L. . 

vq  2  r  u 

(18b) 

A  =  v2  L. . 

VV  11 

(18c) 

The  entropy  production  then  becomes,  in  place  of  Eq.  (9) 

^tot  =  '  ^(t)  '  (19) 


Now,  if  we  consider  the  steady-state  phase  change  problem  to  be  defined 
by  fixing  Ap,  then  the  steady  state  is  characterized  by  minimum  entropy 
production  for  this  constraint  (Refs.  15,  18)  which  requires 


or 


AT 

Ap 


qq 


Substituting  from  Eqs.  (14a),  (14b),  (18a),  and  (18b) 


£-  *1=  2 
T  Ap 


1  + 


2(2-q)(r+l) 

a 


(20a) 


(20b) 


(21a) 


1 8 

Prigogine,  J.,  Thermodynamics  of  Irreversible  Processes.  Interscience 
Publishers,  Division  of  John  Wiley  St  Sons,  Inc.,  New  York,  N.  Y., 

1961. 
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If  Eqs.  (21a)  and  (16b)  are  used  to  eliminate  AT  and  j  ,  respectively,  in 
Eq.  (17b),  it  can  be  shown  that 


Equations  (13a)  and  (13b),  (17a)  and  (17b)  and  (21a)  and  (21b)  are  alternate 
expressions  relating  the  mass  and  energy  fluxes,  allowing  for  both  a  tempera¬ 
ture  and  a  pressure  discontinuity  as  well  as  accounting  for  cross -coupling 
effects  through  the  transport  coefficients  or  A^.  Equation  (21b) 
explicitly  shows  the  relationship  of  the  irreversible  thermodynamics  model 
to  the  Knudsen- Langmuir  model.  Equations  (17a)  and  (17b)  are  in  the  form 
to  be  extended  to  the  case  of  multiple  species  with  chemical  reactions. 

To  employ  the  formulation  developed  in  this  section  in  place  of  the 
Knudsen- Langmuir  model,  one  replaces  Eq.  (la)  with  Eq.  (21b)  and  then 
carries  along  Eq.  (21a)  which  defines  the  additional  unknown  for  this  model, 
i.e.,  the  gas  phase  temperature  T  =  T  +  AT.  Since  C.,  is  the  dominant 

g  w  J 

species  in  the  gas  phase,  the  restrictions  to  C^-only  is  a  reasonable  first 
approximation. 

C.  IRREVERSIBLE  THERMODYNAMICS  MODE  L- -  MU  LT  IP  LE 

SPECIES  WITH  CHEMICAL  REACTIONS 

The  development  of  the  model  for  this  case  is  largely  analogous  to  the 
model  development  outlined  above  for  the  single  component  case.  In  place  of 
Eqs.  (17a)  and  (17b),  which  define  a  2  X  2  matrix  for  j  and  j^,  one  obtains 
in  the  case  of  frozen  gas  phase  chemistry  for  n  species  an  n  X  n  matrix  that 
defines  j  and  (n-1)  independent  species  mass  fluxes  in  terms  of  AT  and 
(n-1)  independent  discontinuities  of  species  chemical  potential  A/i..  For 
the  case  of  n  species  with  gas  phase  chemical  reactions,  one  obtains  a 
(2n-l)  X  (2n-l)  matrix  in  terms  of  the  above  frozen  case  fluxes  and  driving 


F 


potentials  plus  (n-1)  independent  chemical  reaction  rates  and  (n-1) 

chemical  affinities  A..  A  report  describing  this  work  is  being  prepared 

J 

(Ref.  19).  However,  except  for  the  discussion  of  the  effective  heat  of 
ablation,  all  the  results  presented  in  this  paper  were  obtained  either  with 
the  Knudsen- Langmuir  multiple  species  model  or  the  irreversible  thermo¬ 
dynamics  single  component  model. 


T9 

Baker,  R.  L.,  "An  Irreversible  Thermodynamics  Model  for  Graphite 
Sublimation  Including  Multiple  Species  with  Chemical  Reactions,  " 

(In  preparation),  The  Aerospace  Corporation,  El  Segundo,  Calif. 
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III.  RESULTS  AND  DISCUSSION 


A.  KNUDSEN- LANGMUIR  MODEL 

A  complete  description  of  results  calculated  with  this  representation  is 
given  in  Ref.  13.  In  this  section,  a  summary  of  the  results  is  given,  empha¬ 
sizing  the  basis  for  the  need  of  a  more  complete  model. 

The  initial  calculations  performed  were  for  an  external  pressure  of 
one  atmosphere  using  the  1969  JANNAF  thermochemical  property  data 
(Ref.  11)  for  carbon  species  C^-Cj..  Equations  (1)  through  (7)  were  solved 
simultaneously  for  mass  loss  m  and  wall  temperature  T^  as  a  function  of 
the  external  radiation  heat  flux.  The  experimental  data  of  Lundell  and  Dickey 
(Ref.  5)  was  used  to  select  vaporization  coefficients  that  allowed  the  calcu¬ 
lated  results  to  simultaneously  match  both  the  experimental  mass  loss  rate  as 
well  as  the  surface  temperature.  Other  thermochemical  data  bases  (Refs. 

10,  12)  were  then  employed,  and  it  was  found  that  different  values  of  the  species 
vaporization  coefficients  were  required  to  obtain  agreement  of  the  calculated 
results  with  the  experimental  data. 

The  results  of  these  calculations  are  shown  in  Fig.  2  along  with  the 
species  vaporization  coefficient  values  deduced  for  each  thermochemistry 
model.  The  or  values  used  with  the  Dolton,  et  a  I.  ,  thermochemical  data  are 
the  "nominal"  values  selected  in  Ref.  10.  The  unit  vaporization  coefficients 
deduced  for  the  JANNAF  thermochemical  data  are  in  agreement  with  a  prior 
conclusion  of  Lundell  and  Dickey  (Ref.  5).  We  see  from  Fig.  2  that  resolu¬ 
tion  of  a  unique  set  of  thermochemical  data  and  vaporization  coefficients 
cannot  easily  be  accomplished  using  the  Knudsen- Langmuir  model. 

The  major  reason  for  this  nonuniqueness  is  illustrated  in  Fig.  3.  The 
sum  of  the  species  partial  pressures  and  the  sum  of  the  species  vapor 
pressures  are  shown  for  a  wide  range  of  radiation  heat  flux  and  an  ambient 
pressure  of  one  atmosphere.  The  pressure  discontinuity  (£p^  -Xp^) 
increases  dramatically,  moving  from  the  JANNAF  to  the  Palmer  and  then  to 
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the  Dolton  thermo  chemical  data  predictions.  However,  since  the  mass  loss, 
as  given  by  the  Knud  sen- Langmuir  equation  (lb),  is  proportional  to  both  ct^ 
and  (p.V  -  p. ),  smaller  vaporization  coefficient  values  with  large  pressure 
differences  can  compensate  for  unit  vaporization  coefficients  with  small 
pressure  differences. 

The  inability  to  select  a  unique  set  of  thermo  chemical  data  and 
vaporization  coefficients,  while  disturbing  from  a  theoretical  viewpoint, 
appears  from  Fig.  2  to  be  not  too  serious  from  a  practical  viewpoint  since 
all  three  data  sets  reasonably  predict  the  experimental  data.  However,  the 
calculations  performed  have  assumed  that  the  radiation  heat  flux  reaching 
the  surface  is  known.  In  actuality,  the  flux  reaching  the  surface  depends 
upon  radiation  transport  phenomena  occurring  in  the  gas  phase  in  the  path 
of  the  input  radiation  flux.  These  phenomena  in  turn  depend  on  carbon 
species  concentrations  which  are  shown  for  different  thermochemical  models 
in  Fig.  4. 

The  species  mole  fractions  predicted  using  JANNAF  data  when  com¬ 
pared  with  the  Dolton  data  predictions  differ  by  a  factor  of  at  least  two  for 
all  species  in  some  range  of  qrad«  and  the  difference  may  be  as  much  as  a 
factor  of  thirty.  These  uncertainties  provided  the  major  impetus  for 
development  of  the  irreversible  thermodynamics  model.  Comparison  of 
species  concentration  data  with  the  predicted  results  in  Fig.  4  could  poten¬ 
tially  allow  selection  of  a  unique  data  set  using  the  Knud  sen  -  Langmuir 
representation.  However,  as  shown  below,  the  irreversible  thermodynamics 
formulation  provides  additional  useful  information. 

B.  IRREVERSIBLE  THERMODYNAMICS  MODEL 

The  results  discussed,  using  this  model,  in  the  present  report  are 
limited  to  those  in  which  the  carbon  vapor  phase  is  assumed  to  be  homo¬ 
geneous  and  composed  only  of  carbon  species  C^  molecules.  Because  of 
this,  predicted  species  concentration  results  cannot  be  given.  However,  the 
general  nature  of  results  calculated  using  this  extended  model  can  be 
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presented.  In  particular,  new  and  unique  results  obtainable  only  with  the 
irreversible  thermodynamics  model  can  be  illustrated. 


Following  the  procedure  established  with  the  Knudsen- Langmuir  model, 
vaporization  coefficient  values  for  at  were  determined  for  each  of  the 
three  thermo  chemical  data  sets  previously  considered,  by  again  matching 
the  Lundell  and  Dickey  mass  loss  rate-surface  temperature  data.  Appro¬ 
priate  values  of  <*3  were  found  to  be  1.  0,  0.  21  and  0.  062  for  the  JANNAF, 
Palmer  and  Dolton  data,  respectively.  The  predicted  pressure  discontinuities 
Ap  were  found  to  be  similar  to  those  shown  in  Fig.  3.  The  new  variable 
calculated  with  the  present  formulation  is  the  gas  phase  temperature  T  . 

The  predicted  temperature  discontinuities  AT  =  T^  -  T^  are  shown  as  a 
function  of  the  calculated  T^  for  an  ambient  pressure  of  one  atmosphere  in 
Fig.  5.  The  absolute  magnitude  of  AT  is  largest  for  the  JANNAF  data  for 
which  Cl^ssl.O  and  smallest  for  the  Dolton  data  for  which  ct^  -  0.  062  [see 
Eq.  (21a)]. 


Due  to  the  reduction  in  T  below  T  ,  it  is  possible  in  some  cases  for 

g 

the  partial  pressure  of  C3  at  the  ablating  surface  to  be  greater  than  the  vapor 
pressure  of  at  the  gas  phase  temperature,  i.  e,  the  vapor  is  supersaturated. 
Calculated  supersaturation  ratios  are  shown  as  a  function  of  wall  tempera¬ 
ture  and  ambient  pressure  for  the  JANNAF  data  set  in  Fig.  6.  Corresponding 
values  of  qracj  can  be  obtained  from  Fig.  7. 

It  cannot  be  easily  determined  theoretically  if  a  supersaturation  ratio 
of  2-3  is  great  enough  to  cause  condensation  to  solid  carbon.  It  appears  that 
from  the  standpoint  of  the  homogeneous  nucleation  theory,  carbon  species 
super  saturation  levels  much  greater  than  this  are  required  (Ref.  20).  How¬ 
ever,  nucleation  sites  such  as  free  electrons  could  considerably  reduce  the 
required  supersaturation  level.  The  point  is  that  the  present  model  predicts 
that  condensation  may  occur,  and  the  presence  of  solid  carbon  particles  in 


20 

Donn,  B.,  N.  C.  Wickramasinghe,  J.P.  Hudson,  and  T.P,  Stecher,  "On 
the  Formation  of  Graphite  Grains  in  Cool  Stars,  "  Astrophysical  J..  153, 
August  1968,  pp,  451-464. 


Fig.  5.  Predicted  Temperature  Discontinuities, 
Irreversible  Thermodynamics  Model 
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Fig.  6.  Predicted  Supersaturation  Ratios  for 
JANNAF  Chemistry,  a,  =  l.o 


the  flow  would  further  complicate  the  gas  phase  radiation  transport  predic¬ 
tion  problem.  An  additional  important  point  is  that  super  saturation  is  pre¬ 
dicted  to  occur  only  for  the  JANNAF  data  set  which  exhibits  the  largest 

reduction  in  T  below  T  . 

g  w 

C.  GENERAL  RESULTS 

The  above  calculated  results  are  unique  to  the  irreversible  thermo¬ 
dynamics  formulation.  Predicted  mass  loss  as  a  function  of  ambient  pres¬ 
sure  and  qracj  for  the  JANNAF  C^-only  irreversible  thermodynamics  model 

m  0 

is  shown  in  Fig.  7.  For  given  clra(j»  m  is  relatively  insensitive  to  the  ambient 
pressure  level.  This  indicates  that  the  effective  heat  of  ablation  Q*  =  qra{j  /m 
does  not  vary  significantly.  As  pg  decreases,  for  qra{j  constant,  the  flow 
from  the  surface  eventually  becomes  sonic.  This  means  that  further  reduc¬ 
tions  in  pe  have  no  effect.  Therefore,  a  given  point  on  the  solid  line  in 
Fig.  7  at  a  particular  value  of  p^  should  be  interpreted  as  representing  all 
pressures  less  than  p^  also. 

A  given  surface  temperature  can  be  attained  at  many  values  of  p  and 
•  .  e 
qratj.  Generally,  to  reach  a  specified  temperature,  the  highest  qra(j  required 

occurs  when  the  flow  from  the  surface  is  sonic.  As  the  pressure  increases 

and  the  flow  becomes  subsonic  (constant  T^),  the  required  4ra(j  decreases 

significantly.  This  behavior  is  illustrated  more  clearly  in  Fig.  8.  Figure  8 

shows  that  if  the  melt  temperature  is  low,  i.e.,  3800°K  (Ref.  21),  then 

the  qra<J  required  for  melting  to  occur  may  be  nearly  an  order  of  magnitude 

below  that  required  to  reach  the  more  commonly  accepted  melt  temperature 

of  4200°K  (Ref.  22). 
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Whittaker,  A.G.,  and  P.  L,  Kintner,  "Laser-Heating  Studies  of  Carbon 
Melting  and  Vaporization,  "  American  Ceramic  Society  27th  Pacific  Coast 
Regional  Meeting,  North  Hollywood,  Calif.,  Oct.  23-26,  1974. 

22 

Diaconis,  N.  W.,  E.  R.  Stove,  J.  Hook,  andG.J,  Catalano,  "Graphite 
Melting  Behavior,  "  Re-Entry  and  Environmental  Systems  Division, 

General  Electric  Company,  AFML-TR-71- 119,  July  1971. 


Fig.  8.  Relationship  of  Tw  to  qra(j  at  Several  Pressures. 
Irreversible  Thermodynamics  Model 
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The  question  of  melt  temperature  remains  open  at  the  present  time 
(Refs.  23,  24).  If  melting  does  not  occur,  the  variation  in  the  effective  heat 
of  ablation  Q  ,  predicted  by  all  of  the  various  models  considered  in  this 
report,  is  relatively  small.  This  is  seen  in  Table  1  which  includes  results 
predicted  by  the  5- species  irreversible  thermodynamics  model.  If  melting 
does  occur,  the  Q  value  could  be  decreased  an  order  of  magnitude  by 
mechanical  removal  of  liquid  carbon  (Ref.  21). 


Table  1.  Q*  Values 


Model 

& 

Range  of  Q 
(MJ/kg) 

Knudsen-  Langmuir /JANNAF 

31.  0  -  33.  1 

Knudsen- Langmuir /Palmer 

28.4  -  31.3 

Knudsen-  Langmuir /Dolton 

26.  9  -  32.3 

Present-C3  Only/JANNAF 

27.4  -  28.  1 

Present-C^  Only/Palmer 

27.4  -  28.4 

Present-C^  Only /Dolton 

29.3  -  30.3 

Present-5  Species/JANNAF 

30.  5  -  32.5 
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High  Temperature  Science.  8,  1976,  pp.  81-97. 
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Whittaker,  A.  Greenville,  and  P,  K.  Kinter,  "Comments  on  the  Solid- 
Liquid- Vapor  Triple  Point  of  Carbon,  "  High  Temperature  Science. 

9(1),  March  1977.  - 
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The  results  presented  herein  are  preliminary  in  many  respects.  It  is 
anticipated  that  further  comparisons  of  the  complete  5-species  irreversible 
thermodynamics  model  predictions  with  experimental  data  will  help  facilitate 
the  selection  of  a  unique  set  of  thermochemical  property  data  and  the  correct 
melt  temperature.  The  following  experimental  data  would  aid  considerably  in 
this  work:  (1)  carbon  species  concentration  measurements  above  radiation 
heated  graphite  surfaces,  (2)  measurements  of  the  reduction  in  the  gas  phase 
temperature  below  the  surface  temperature,  (3)  detection  of  condensation  to 
solid  carbon  particles,  and  (4)  detection  of  surface  melting. 
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IV.  SUMMARY  AND  CONCLUSIONS 


It  has  been  demonstrated  by  comparing  predicted  results  with  experi¬ 
mental  data  that  a  unique  set  of  vaporization  coefficients  and  species  ther- 
moehemical  data  cannot  be  readily  selected  using  a  Knudsen- Langmuir  type 
multiple  species  model.  Because  of  this,  a  large  uncertainty  presently 
remains  regarding  the  concentration  of  carbon  species  C^-C^  above  a  carbon 
surface  exposed  to  a  high  radiation  heat  flux.  This  uncertainty  will  be  very 
important  to  consider  if  gas  phase  radiation  transport  predictions  are 
required. 

In  an  attempt  to  help  resolve  this  uncertainty,  an  irreversible  thermo¬ 
dynamics  phase  change  model  has  been  adapted  to  this  problem  in  order  to 
incorporate  effects  neglected  by  the  Knudsen- Langmuir  representation. 

Preliminary  results  calculated  with  the  new  model,  assuming  gaseous 
carbon  to  be  composed  of  molecules  only,  indicate  significant  differ¬ 
ences  in  the  predicted  temperature  discontinuity  between  the  subliming  sur¬ 
face  and  the  gas  phase,  depending  on  the  thermochemical  data  used.  The 
largest  temperature  differences  are  predicted  for  the  JANNAF  data,  and 
the  reduced  gas  temperatures  in  this  case  result  in  predicted  gas  phase 
supersaturation  ratios  of  as  much  as  2-3.  For  the  other  thermochemical 
data  considered,  no  supersaturation  is  predicted. 

Assuming  that  the  graphite  does  not  melt  with  subsequent  liquid  removal, 
the  predicted  effective  heat  of  ablation  Q  is  30  ±  3  MJ/kg  for  all  the  dif¬ 
ferent  model  types  and  thermochemical  data  sets  used.  However,  if  melting 
with  liquid  removal  occurs,  this  could  be  reduced  considerably.  The  pre¬ 
dicted  radiation  flux  required  to  reach  various  postulated  melt  temperatures 
is  given  as  a  function  of  ambient  pressure  in  Fig.  8. 

Additional  work  employing  the  irreversible  thermodynamics  phase 
change  model  considering  multiple  species  with  chemical  reactions  is  in 
progress. 
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